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We give a short update of our research program on nonequilibrium statistical field theory 
fT") . applied to quantum processes in the early universe and black holes, as well as the develop- 

' ment of stochastic gravity theory as an extension of semiclassical gravity and an intermediary 

, in the 'bottom-up' approach to quantum gravity. 
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O*. PROLOGUE 

<: 

\Q ■ (BLH) "It is no secret that the message of this conference is 'NOISE IS GOOD'. In this talk I want to show 

that not only is noise good, it is absolutely essential... " . These words were not uttered in this conference, as 
it would have been acausal, but ten years ago, at a workshop devoted to the nascent yet fascinating subject 
of Fluctuations and Order [1]. After a decade, the tenor and objectives of this statement are even closer to 
reality, as witnessed by the strong focus and great variety of this symposium. 
J^-j I substantiated this claim then by an enumeration of the many processes in gravitation and cosmology 

, where quantum noise and fluctuations play an active role: 

1. Particle creation as parametric amplification of vacuum fluctuations 

2. Thermal radiance from accelerated observers and black holes as fluctuation-dissipation phenomena 

3. Entropy generation from quantum stochastic and kinetic processes 

4. Phase transitions in the early universe as noise-induced processes 
O 5. Galaxy formation from primordial quantum fluctuations 

6. Anisotropy dissipation from particle creation as backreaction processes 
i—t , 7. Dissipation in quantum cosmology and the issue of the initial state 

8. Decoherence, backreaction and the semiclassical limit of quantum gravity 

9. Stochastic spacetime and continuum limit, gravity as an effective theory 

10. Topology change in spacetime and loss of quantum coherence problems 

11. Gravitational entropy, singularity and time asymmetry 

12. 'Birth' of the universe as a spacetime fluctuation and tunneling phenomenon 

The above list was prepared for cosmological issues. Further developments in the last ten years focusing 
on the effects of noise and fluctuations using the statistical field theory approach we have developed (for 
a review, see [2,3]) include, related to Topic 1 : Preheating in post-inflationary cosmology [4,5]; Topic 2: 
Thermal and near-thermal radiance in detectors, moving mirrors, black holes and cosmology [6-9]; Topic 3: 
Correlation entropy [10,11]; Topic 4: Defect formation [12], tunneling induced by quantum and thermal noise 
fluctuations [13,14]. Topic 5 is presently under pursuit [16] and partially summarized in the last section of this 
paper. Topic 9 includes preliminary work on wave propagation in stochastic spacetimes [17] and mesoscopic 
fluctuations [18]. 

In black hole physics, focusing again on fluctuations, one can mention topics on black hole fluctuations and 
backreaction [19,20] (see references therein for other related work) and the energetics and dynamics of black 
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hole phase transition [21-23]. The wish list would include: applications of stochastic gravity to the statistical 
mechanical definition of entropy and statistical field theory of black hole noncquilibrium thermodynamics. 

Using statistical quantum field theory , Topics 1, 2, 3, 6 are quite well understood, 7, 8 only partially, and 
Topics 4, 5, 9 are currently under pursuit. Certain concepts in Topics 10-12 may not even be well-defined 
(e.g., what does the 'birth' of the universe mean?), but the newly established theory of stochastic gravity (for 
a review, see [24-27], an ongoing development, see [28]) may offer alternative ways to address these issues, 
as well as provide an intermediary towards quantum gravity (our definition is different from the ordinary, 
see [29]). 

In this talk I will focus on the theory of noise in quantum fields and how quantum fluctuations could 
have played an active and even decisive role in many fundamental processes in cosmology and gravitation, 
especially near the Planck time (10~ 43 sec from the Big Bang). I will describe how stochastic gravity 
theory can be understood easily from an open system conceptual framework, and give two examples of its 
applications: fluctuations and backreaction in black holes, and structure formation in the early universe. 

The Planck time is the time when many familiar features of spacetime depicted by Einstein's theory of 
general relativity give way to an as-yet-unknown quantum theory of gravity depicting the microstructure 
of spacetime. (Some theoreticians believe the superstring theory is the answer.) Just below the Planck 
energy we believe the universe can be adequately described by a semiclassical theory of gravity [30-32], 
where quantized matter fields coexist with a classical spacetime. Many qualitative changes are believed to 
have taken place at this energy scale, amongst them the formation of spacetime depictable as a manifold, 
the emergence of time, the creation of particle pairs from the vacuum, the growth of fluctuations as seeds 
for galaxies, and possible phase transitions and the ensuing entropy generation processes. It is also the 
cross-over point of quantum to classical and micro to macroscopic transitions. 1 

In Section I, we first explain the origin and nature of noise in quantum systems interacting with an 
environment, using the influence functional method. (For non-Ohmic bath at low temperatures, colored 
noise and nonlocal dissipation would appear; and for nonlinear coupling multiplicative noise is generally 
expected. A generalized fluctuation-dissipation relation for these systems can be proven, and the stochastic 
(master, Langevin and Fokker-Planck) equations derived, depicting the dissipative dynamics of the open 
system under the influence of noise. We then discuss stochastic gravity theory in Sec. 2 and the two 
examples in Sections 3 and 4. 

I. QUANTUM FLUCTUATIONS IN OPEN SYSTEMS 

We begin by describing how the concepts of quantum open systems can be of use in the treatment of 
statistical mechanical problems involving quantum fields. We start with two subsystems A and B. When the 
precise information of subsystem B is not required, but only its averaged effect on subsystem A is of interest, 
one can coarse-grain B and include its averaged effect on A, which involves finding the backreaction on A. In 
so doing A is rendered an open system, with B acting as its environment. For the analysis of open systems 
with backreaction from the environments the influence functional (IF) formalism of Feynman and Vernon 
[35,36] proves useful. Let us first use a simple example from quantum mechanics to illustrate the idea and 
the method. We then show how we can use this framework to address issues in gravitation and cosmology. 



1 lt is for this reason that I think one can view general relativity as geometro-hydrodynamics [34], semiclassical gravity 
as a mesoscopic physics [33], and take a kinetic theory approach to quantum gravity [29]. The only difference is that 
instead of dealing with the quantum to classical and micro- to macro- transition in the state of matter and fields we 
are dealing with the corresponding issues for spacetime and geometry. 
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A. Stochastic Effective Action in Quantum Open Systems 



We begin with a brief schematic summary of the IF formalism as applied to a simple system. Consider 
a system S, described by the degrees of freedom x, interacting with an environment E, described by the 
degrees of freedom q. 2 The full closed quantum system S + E is described by a density matrix p(x, q; x' , q' , t). 
If we are interested only in the state of the system as influenced by the overall effect, but not the precise 
state of the environment, i.e, the dynamics of the open system, then the reduced density matrix p r (x, x',t) = 
J dq p(x, q; x' , q, t) would provide the relevant information. (The subscript r stands for reduced.) Assuming 
that the action of the coupled system decomposes as S = S s [x] + S e [q] + Si n t [x, q] , and that the initial density 
matrix factorizes (i.e., takes the tensor product form), p(x,q;x',q',ti) = p s (x,x' ,ti)p e (q,q' , t,), the reduced 
density matrix is given by 

+00+00 , 

Pr (x,x',t)= J dx, J dx\ Dx P Dx' e i(Ss[x]-S s [x']+S IF [ X ,x',t]) ^fa^. ) ( L 1) 

— 00 — 00 

where Sif is the influence action related to the influence functional T defined by 

/ris ris 
dqfdqidq'i Dq Dq' e ^SM+S mt [x., q ]-S c [ q ]-S int \x ,g ]) p ^ qu ^ ( L2 ) 

Sif in general is complex. Retaining only quadratic terms (an approximation which covers many of the 
interesting applications that we will consider later), we may write 



S IF (x,x') = [ dtdt' {\{x - x'){t)D{t,t'){x + x'){t') + Ux - x')(t)N{t,t'){x - x'){t')} 
J Z 2 



(1.3) 



where D and N stand for the real dissipation and noise kernels respectively. Note that in this quadratic 
order approximation, the influence action Sif{x, x') is related to the closed-time-path(CTP) or in-in effective 
action (for details on the CTP effective action see [38]) Tctp[x, x'] through 

T C tp[x,x'} = S[x] - S[x'] + S IF [x,x'\. (1.4) 

The equation of motion obtained from the CTP effective action for the expectation values is clearly seen to 
be real and causal [38] . It reads 



r C Tp[x,x'] 



5x(t) 



(1.5) 



x'—x—x 



From the influence functional a Langevin equation for the system dynamics may be derived by a formal 
procedure, first introduced by Feynman and Vernon [35], which consists of introducing a Gaussian stochastic 
source £(t) with (£(£)){ = and ( = -^(M') an d defining an improved or stochastic effective action 

as 

S eff [x, x'-^} = S s [x] - S s \x'] + TZS IF [x, x'} + £(x - x') (1.6) 
such that (^e lSe f/[ x ' x ''^ — e ,rcTf, ^ I 'l This leads to equations of motion with a stochastic force: 



dS ef f[x,x';^] 



5x 



n ■ i 4.1 <5 r CTp[x,x';^] 

or, equivalently, 



(5x 



0. (1.7) 



2 We are labeling the degrees of freedom of the system and the environment by single letters x and q with the 
understanding that they can represent multiple or even infinite degrees of freedom, e.g. corresponding to a field [37]. 
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The equation of motion obtained from (1.7) using (1.6) is 



dS s 
dx(t) 



+ 



(1.8) 



where D(t,t') = —d t 'j(t,t'). Being now in the form of a Langevin equation, the physical meaning of the 7 
and N kernels in Eq. (1.8) becomes clearer. Both the terms involving 7 and £ represent the backreaction 
of the environment on the system. However, 7 (or more properly the odd part of 7) is associated with 
dissipation and £ is a stochastic noise term associated with random fluctuations of the system, exactly as 
the terms are interpreted in the context of Brownian motion. Averaging (1.8) over the noise using the 
appropriate probability distribution will give the semiclassical equation of motion for the mean value of x. 
Since the noise and dissipation arise by considering a subsystem within a closed system (as is done here, as 
opposed to being put in by hand), they are in general related by a set of generalized fluctuation-dissipation 
relations (FDR), which can be represented by a linear, non-local relation of the following form, provided that 
the Hamiltonian for the environment and the system-environment interaction are time independent and the 
initial state of the environment is stationary with respect to the Hamiltonian of the environment: 



To keep the discussion simple, we have written the noise and dissipation kernels in terms of single scalar 
functions. However, the method is general enough to encompass multiple noise and dissipation kernels and 
cases where the kernels are tensorial, as in the stochastic gravity theory discussed later. 



Stochastic semiclassical gravity [39] of the 90's is a theory naturally evolved from semiclassical gravity [32] 
of the 80's and quantum field theory in curved spacetimes [30,31] of the 70's. Whereas semiclassical gravity is 
based on the semiclassical Einstein equation with sources given by the expectation value of the stress-energy 
tensor of quantum fields, stochastic semiclassical gravity is based on the Einstein-Langevin equation, which 
has in addition stochastic sources with correlation functions characterized by the noise kernel, the noise 
kernel. The noise kernel is the vacuum expectation value of the (operator-valued) stress-energy bi-tensor 
which describes the fluctuations of quantum matter fields in curved spacetimes. 



The first stage in the road to stochastic gravity begins with quantum field theory in curved spacetime, which 
describes the behavior of quantum matter fields propagating in a specified (not dynamically determined by 
the quantum matter field as a source) background gravitational field. For a scalar field <f> it obeys the 
wave equation (□ + m 2 )(f>{x) = where □ is the Laplace-Bcltrami operator, which contains the imprint of 
the curvature of the background spacetime. In this theory the gravitational field is given by the classical 
spacetime metric determined from classical sources by the classical Einstein equations, and the quantum 
fields propagate as test fields in such a spacetime. For time dependent spacetime geometry it may not be 
possible to define a physically meaningful vacuum state for the quantum field at all times. Assuming that 
one defines a vacuum state at some initial time, the vacuum state at a latter time will differ from that 
defined initially because particles are created in the intervening time. An important process described by 
quantum field theory in curved spacetime is indeed particle creation from the vacuum (and effects of vacuum 
fluctuations and polarizations) in the early universe [41] and Hawking radiation in black holes [41]. 

The second stage in the description of the interaction of gravity with quantum fields is back-reaction, i.e., 
the effect of quantum fields on the spacetime geometry. The dynamic classical spacetime metric creates 
particles of the quantum field and these in turn provide a backreaction on the spacetime metric which alters 
its dynamics in response. One assumes a general class of spacetime where the quantum fields live in and act 




(1.9) 



II. STOCHASTIC GRAVITY AND METRIC FLUCTUATIONS 



A. From Semiclassical to Stochastic Gravity 
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on, and seek a solution which satisfies simultaneously the Einstein equation for the spacetime and the field 
equations for the quantum fields. The Einstein equation which has the expectation value of the stress-energy 
operator of the quantum matter field as the source is known as the semiclassical Einstein equation: 

G» v {g a p) = K{% v ) q (III) 

where 7j,„ is the stress-energy tensor operator of, say, a free scalar field <fi, G^ u is the Einstein tensor, 
k = 8ttGn and Gn is Newton's constant. Here ( ) q denotes the expectation value taken with respect to some 
quantum state compatible with the symmetries of the background spacetime, a classical object. The theory 
obtained from a self-consistent solution of the geometry of the spacetime and the quantum field is known as 
semiclassical gravity. Incorporating the backreaction of the quantum matter field on the spacetime is thus 
the central task in semiclassical gravity. 

Studies of the semiclassical Einstein equation for the backreaction problems have been carried out in 
the last two decades by many authors for cosmological and black hole spacetimes. A well-known example 
of semiclassical gravity is the damping of anisotropy in Bianchi universes (which is the basis of chaotic 
cosmology in the 70's) by the backreaction of vacuum particle creation, and inflationary cosmology [42-44] 
of the 80's driven by a constant vacuum energy density source such as the expectation value of a Higgs field. 

In analogy with the open system dynamics described in Section I A, Eq.(II.l) is equivalent to Eq.(1.5) 
where the degrees of freedom x of the system are identified with the metric g a p and those of the environment 
q are identified with the scalar field 4>(x). However, from the discussion in the last section it is also clear 
that Eq. (1.5) , and hence also the semiclassical Einstein Eq. (II. 1) results on averaging the full Langevin- 
type Eq. (1.7) over noise. Thus the semiclassical Einstein equation incorporates the dissipation but misses 
out the fluctuation aspect of the backreaction. The recognition of this crucial point [45] ushered in a new 
theory known as stochastic semiclassical gravity, (or in short, stochastic gravity, as there is no confusion 
in this context as to where the stochasticity originates). Aided by the concept of open systems and the 
techniques of the influence functional and the Closed Time Path (CTP) effective action, stochastic gravity is 
the new framework for the consideration of backreaction because it encompasses fluctuations and dissipation 
(from particle creation and other quantum field processes) on the same footing. Spacetime dynamics is 
now governed by a stochastic generalization of the semiclassical Einstein equation known as the Einstcin- 
Langevin equation, the analog of Eq. (1.8) in the context of semiclassical gravity (SCG). Schematically the 
Einstein-Langevin equation takes on the form 

G i1v {x)=k(T; v +T*1), T^^(f^) q + T^ (11.2) 

Here, T° is due to classical matter or fields, {T^ v ) q is the expectation value of the stress tensor of the 
quantum field, and T/*£ is a new stochastic term which is related to the quantum fluctuations of T^ v for the 
state of the field under consideration. Taking the average of (II. 2) with respect to the noise distribution will 
lead to the conventional semiclassical Einstein equation. It is in this context that SCG is regarded as a mean 
field theory 

The fundamentals of this new theory were developed via two approaches: the axiomatic and the functional. 
The axiomatic approach is useful to see the structure of the theory from the framework of semiclassical 
gravity, showing the link from the mean value of the stress-energy tensor to its correlation functions. The 
functional approach uses the Feynman- Vernon influence functional and the Schwinger-Keldysh closed-timc- 
path effective action methods which are convenient for computations. It also brings out the open systems 
concepts and the statistical and stochastic contents of the theory such as dissipation, fluctuations, noise and 
decoherence. There was also theoretical work on the properties of the stress energy bi-tensor and its vacuum 
expectation value, the noise kernel. See, e.g., [46-52]. For a broader exploration of ideas and issues based 
on this theory read the reviews [24,27]; For a pedagogical introduction with applications, see [26,25]. 

Thus with the aid of the open system viewpoint it is easy to see that stochastic gravity is a natural 
extension of the well-established semiclassical gravity theory and a useful framework for the considerations 
of fluctuations in quantum matter fields and dissipative dynamics of classical spacetimes, including metric 
fluctuations. Stochastic gravity can address many important issues related to nonequilibrium quantum field 
processes in curved spacetimes and find applications to many problems in gravitation and cosmology. 
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B. Stochastic Gravity in relation to Quantum Gravity 

Before embarking on the discussion of some applications of stochastic gravity, let us illustrate the theory 
with a simple toy model which minimizes the technical complications. The model will be useful to clarify 
the role of the noise kernel and illustrate the relationship between the semiclassical, stochastic and quantum 
descriptions. Let us assume that the gravitational equations are described by a linear field h(x) whose source 
is a massless scalar field <j)(x) which satisfies the Klein-Gordon equation in flat spacetime U(p(x) = 0. The 
field stress-energy tensor is quadratic in the field, and independent of h(x). The classical gravitational field 
equations will be given by 

Uh{x) = kT(x), (II.3) 

where T(x) is the (scalar) trace of the stress-energy tensor, T(x) — d a <fi(x)d a <j)(x) and k = 16itG, where 
G is Newton's constant. Note that this is not a self-consistent theory since 4>(x) does not react to the 
gravitational field h(x). We should also emphasize that this model is not the standard linearized theory of 
gravity in which T is also linear in h(x). It captures, however, some of the key features of linearized gravity. 
In the Heisenberg representation the quantum field h(x) satisfies 

Uh{x) = kT(x). (II.4) 

Since T(x) is quadratic in the field operator <j>(x) some regularization procedure has to be assumed in order 
for (II. 4) to make sense. Since we work in flat spacetime we may simply use a normal ordering prescription 
to regularize the operator T(x). The solutions of this equation, i.e. the field operator at the point x, h(x), 
may be written in terms of the retarded propagator G(x, y) as, 

h(x) = {x)+K J dx'G(x, x')f(x'), {11.5) 

where (x) is the free field which carries information on the initial conditions and the state of the field. 
From this solution we may compute, for instance, the symmetric two point quantum correlation function 
(the anticommutator) 

i({h(x),h(y)}) = l({hW(x),h^(y)}) + ^- J J dx'dy'G(x,x')G(y,y'){{f(x'),f(y')}), (II.6) 

where the expectation value is taken with respect to the quantum state in which both fields 4>(x) and h{x) 
are quantized. (We assume for the free field, (h^) = 0.) 

We can now consider the semiclassical theory for this problem. If we assume that h(x) is classical and 
the matter field is quantum the semiclassical theory may just be described by substituting into the classical 
equation (II. 3) the stress-energy trace by the expectation value of the stress-energy trace operator (T(x)), 
in some quantum state of the field <f>(x). Since in our model T(x) is independent of h(x) we may simply 
renormalize its expectation value using normal ordering, then for the vacuum state of the field (f>(x), we 
would simply have (T(x))o = 0. The semiclassical theory thus reduces to 

Uh{x) = k(T(x)). (11.7) 

The two point function h(x)h(y) that one may derive from this equation depends on the two point function 
{T(x))(T(y)) and clearly cannot reproduce the quantum result (II. 6) which depends on the expectation value 



3 In this article we use the (+,+,+) sign conventions of Ref. [53], and units in which c = Ti — 1. 
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of the two point operator ({f(x), T ($/)}). That is, the semiclassical theory entirely misses the fluctuations 
of the stress-energy operator T(x) . 

Let us now see how we can extend the semiclassical theory in order to account for such fluctuations. 
The first step is to characterize these fluctuations. For this, we introduce the noise kernel as the physical 
observable that measures the fluctuations of the stress-energy operator T. Define 

N{x,y) = l -{{i{x),t{y)}) (H.8) 

where t(x) — T[x) — (T(x)). The bi-scalar N(x, y) is real and positive-scmidcfinite, a consequence of t being 
self-adjoint. A simple proof can be given as follows. Let \ip) be a given quantum state and let Q be a 
self-adjoint operator, — Q, then one can write (ip\QQ\ip) — (ip\Q^Q\ijj) — \Q\ip)\ 2 > 0. Now let t(x) be 
a self-adjoint operator, then if we define Q — J dxf(x)t(x) for an arbitrary well behaved function f(x), the 
previous inequality can be written as J dxdyf (x)(ip\i(x)i(y)\ip} f (y) > 0, which is the condition for the noise 
kernel to be positive semi-definite. Note that when considering the inverse kernel N" 1 (x,y), it is implicitly 
assumed that one is working in the subspace obtained from the eigenvectors which have strictly positive 
eigenvalues when the noise kernel is diagonalized. 

By the positive semi-definite property of the noise kernel N(x,y) it is possible to introduce a Gaussian 
stochastic field as follows: 

=0, (£(x)t(y)) a =N(x,y). (II.9) 

where the subscript s means a statistical average. These equations entirely define the stochastic process £(x) 
since we have assumed that it is Gaussian. Of course, higher correlations could also be introduced but we 
just try to capture the fluctuations to lowest order. 

The extension of the semiclassical equation may be simply performed by adding to the right-hand side 
of the semiclassical equation (II. 7) this stochastic source £(x) which accounts for the fluctuations of T as 
follows, 

Uh{x) = k (<T (a;)) + £(a;)) . (11.10) 

This equation is in the form of a Langevin equation: the field h(x) is classical but stochastic and the 
observables we may obtain from it are correlation functions for h(x). In fact, the solution of this equation 
may be written in terms of the retarded propagator as, 

h(x) = (x)+k J dx'G(x, x 1 ) ((TV)) + . 

from where the two point correlation function for the classical field h(x), after using the definition of £(x) 
and that (h^°\x)) s — 0, is given by 

(h(x)h(y)) s = (h(°\x)hW(y)) a + y // dx'dy'G(x,x')G(y,y')({f(x'),f(y')}). (11.12) 

Note that in writing (. . .) s here we are assuming a double stochastic average, one is related to the stochastic 
field £(x) and the other is related to the free field h^(x) which is assumed also to be stochastic with a 
distribution function to be specified. 

Comparing (II. 6) with (11.12) we see that the respective second term on the right-hand side are identical 
provided the expectation values are computed in the same quantum state for the field 4>{x) (recall that we 
have assumed T(x) does not depend on h(x)). The fact that the field h(x) is also quantized in (II. 6) does 
not change the previous statement. The nature of the first term on the right-hand side of Eqs. (II. 6) and 
(11.12) is different: in the first case it is the two point quantum expectation value of the free quantum field 
M°) whereas in the second case it is the stochastic average of the two point classical homogeneous field h^°\ 
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which depends on the initial conditions. Now we can still make these terms equal to each other if we assume 
for the homogeneous field a Gaussian distribution of initial conditions such that 

(ft<°> (*)/»<%)>. = \{{hM(x)rt 0) (v)})- (H.13) 

This Gaussian stochastic field (x) can always be defined due to the positivity of the anti-commutator. 
Thus, under this assumption on the initial conditions for the field h(x) the two point correlation function 
of (11.12) equals the quantum expectation value of (II. 6) exactly. An interesting feature of the stochastic 
description is that the quantum anticommutator of (II. 6) can be written as the right-hand side of equation 
(11.12), where the first term contains all the information on the initial conditions for the stochastic field h(x) 
and the second term codifies all the information on the quantum correlations of the source. This separation is 
also seen in the description of some quantum Brownian motion models which are typically used as paradigms 
of open quantum systems [13,14]. 

It is interesting to note that in the standard linearized theory of gravity T(x) depends also on h(x), both 
explicitly and also implicitly through the coupling of 4>{x) with h{x). The equations are not so simple but it 
is still true that the corresponding Langevin equation leads to the correct symmetrized two point quantum 
correlations for the metric perturbations [48,28]. Thus in a linear theory as in the model just described 
one may just use the statistical description given by (11.10) to compute the symmetric quantum two point 
function of equation (II. 5). This does not mean that we can recover all quantum correlation functions with 
the stochastic description, see Ref. [13] for a general discussion about this point. Note that, for instance, the 
commutator of the classical stochastic field h(x) is obviously zero, but the commutator of the quantum field 
h(x) is not zero for timelike separated points; this is the prize we pay for the introduction of the classical field 
£(x) to describe the quantum fluctuations. Furthermore, the statistical description is not able to account for 
the graviton-graviton effects which go beyond the linear approximation in h{x). 



III. BLACK HOLE FLUCTUATIONS AND BACKREACTION 

As the first example we consider the backreaction of Hawking radiation on black holes [54] with fluctuations, 
i.e., how a quantum field and its fluctuations influence the behavior of the background spacetime. We will 
only sketch the strategy of this program based on stochastic gravity, as detailed calculations are still in 
progress. The formalism described in Section I A will be useful. Here we study the simpler class of problems 
of a quasi-static black hole in quasi-equilibrium (a box is required) with its Hawking radiation described by 
a scalar field. The goal is to obtain a stochastic influence action analogous to (1.6) for this model of a black 
hole coupled to a scalar field. From it one can derive an Einstein-Langevin equation analogous to (1.8). 

We consider the simplest model of this class described by a perturbed Schwarzschild metric, used by York 
[55] to analyze black hole backreaction. We focus on the new aspects of noise and fluctuations, their origin 
and attributes. 

In this model the black hole spacetime is described by a spherically symmetric static nonvacuous metric 
with line element of the following general form written in advanced time Eddington-Finkelstcin coordinates 

ds 2 = g^dx^dx" = -e 2 ^ (l - ^ dv 2 + 2e 2qp dvdr + r 2 dtl 2 (III.l) 

where ip = ip(r) and m = m(r) , v = t + r + 2Mln (^j — l) and dfi 2 is the line element on the two sphere. 
Hawking radiation is described by a massless, conformally coupled quantum scalar field 4> with the classical 
action 

S m [</>, =~\j d n x^—g[g^d^d v 4> + an)R(b 2 } (HI.2) 

where £ (n) = (n is the dimension of spacetime) and R is the curvature scalar of the spacetime it lives 

in. 
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Let us consider linear perturbations = g^ u — gffi off a background Schwarzschild metric g$ with line 
element 

(ds 2 ) = (l - dv 2 + 2dvdr + rW (III.3) 

We look for this class of perturbed metrics in the form given by (III.l), (thus restricting our consideration 
only to spherically symmetric perturbations): 

~ 1 + ep{r), m ~ M[l + ep(r)\ (III.4) 

where = ^aT^;a = A = 90(8 4 )7r 2 . T h is the Hawking temperature. (This particular parametriza- 
tion of the perturbation is chosen following York's [55] notation.) Thus the only non-zero components of /i M „ 
are 

h vv = -({!- —)2ep(r) + 2Me ^ r A h vr = ep(r) (111.5) 
\ r r J 

So this represents a metric with small static and radial perturbations about a Schwarzschild black hole. The 
initial quantum state of the scalar field is taken to be the Hartle-Hawking vacuum, which is essentially a 
thermal state at the Hawking temperature as far as static observers are concerned 4 and it represents a black 
hole in (unstable) thermal equilibrium with its own Hawking radiation. 

The metric perturbation expansion induces a decomposition of the Einstein tensor G M „ ~ G|S + 5G^ V 

where G$ is the Einstein tensor for the background spacetime. The zeroth order solution gives a background 
metric in empty space, i.e, the Schwarzschild metric. 8G^ V is the linear correction to the Einstein tensor in 
the perturbed metric. The semiclassical Einstein equation in this approximation therefore reduces to 

6G^°\h) = K(T^) 

York solved this equation to first order by using the expectation value of the energy momentum tensor 
for a conformally coupled scalar field in the Hartle-Hawking vacuum in the unperturbed (Schwarzschild) 
spacetime on the right hand side and <5G M „ on the left hand side is calculated using (III. 5). This yields 
the corrections to the background metric induced by the backreaction encoded in the functions p(r) and 
p(r), which amounts to the noise-averaged backreaction effects. We are interested in the fluctuations and its 
effects. 

We now derive the CTP effective action for this model, following the treatment of Ref . [56] . Using the 
metric (III. 3) (and neglecting the surface terms that appear in an integration by parts) we have the action 
for the scalar field written perturbatively as 

S m [4>, V] = \J d n x^*) 4, [□ + + + . . .] 0, (IH.6) 

where the first and second order perturbative operators and are given by 

V {2) = --j^{[dp(y/^h^{x))]d v + h^{x)dv^ (III.7) 



4 The Hartle-Hawking vacuum is a pure state which is perceived as vacuum by free-falling observers crossing either 
the black hole or the white hole horizons. However, if one traces out the modes localized in the second asymptotically 
flat region, one ends up with an incoherent density matrix perceived as a thermal state by static observers in the first 
asymptotically flat region. 
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In the above expressions, R^ is the fc-order term in the pertubation h^ v (x) of the scalar curvature R and 
h^v and denote a linear and a quadratic combination of the perturbation, respectively, 

V = K» - \h9$\ V = K h ™ - \hh^ + ±h 2 g$ - \h a ph a Pg$. (III.8) 

From quantum field theory in curved spacetime considerations we take the following action for the gravita- 
tional field (see [47,56] for more details) 

m^{x)R^{x) (\ - 360(£(n) - i) 2 ) R(x)R{x) . 

(III.9) 

The first term is the classical Einstein-Hilbert action and the second term is the countcrtcrm in four dimen- 
sions used to renormalize the divergent effective action. In this action £ 2 p = 16nG, a = (28807T 2 ) -1 and p, is 
the renormalization mass scale. 

We are interested in computing the CTP effective action for the model given by the form (III. 6) for the 
matter action and when the field <j> is initially in the Hartlc- Hawking (HH) vacuum. Since the initial state 
of the field is described by a thermal density matrix at the HH temperature Th 5 , the finite temperature 
CTP effective action (T = 1//3) for this model is given by (for details see [56]) 

^ctpHA = SW[h+,] Sf^[h- V ] ^Tr{\nG ab [h%]}, (111.10) 

where ± denote the forward and backward time path of the CTP formalism and G^ b [h^ v ] is the complete 

2x2 matrix propagator (a and b take ± values: G ++ ,G^ and G__ correspond to the Feynman, Wightman 

and Schwingcr Green's functions respectively) with thermal boundary conditions for the differential operator 

y/-g(°)(D + + H ). The actual form of <5f fc cannot be explicitly given. However, it is easy to 

obtain a perturbative expansion in terms of V^ b \ the fc-order matrix version of the complete differential 
operator defined by V±± = ±V^ and V±^ = 0, and G^ b , the thermal matrix propagator for a massless 
scalar field in Schwarzschild spacetime . To second order G^ b reads 

Gl = (& ~ GlV$G% - GlV$G% + GlV$Giv}?G% + ■■■ (III.ll) 

Expanding the logarithm and dropping one term independent of the perturbation h^ v {x), the CTP effective 
action may be perturbatively written as 

^ctpH v ] = ] - Sf[h- V ] + l -Tr[vl r) G f i + V™Gt_ + if } G^ + - V™Gt_] 

-iTriV^G^V^G^ + V V>GfL_vV>Gt_ - 2V™ G^V™ Gt + ]. (111.12) 

However, unlike the case of Ref. [56] where represented a perturbation about flat space and hence one 
had knowledge of exact "unperturbed" thermal propagators, in this case, since the perturbation is about 
Schwarzschild spacetime, exact expressions for the corresponding unperturbed propagators G^ b [h^ u ] arc not 
known. Therefore apart from the approximation of computing the CTP effective action to certain order 



i p — n — 4 p 

sf v) [g,A = J d»x j-gR{x) + J d n x 



5 As mentioned earlier, this is true for the basis of modes associated to static observers and provided that one is 
not concerned about quantum correlations (entanglement) between the two asymptotically flat regions, so that the 
modes localized in the second region can be traced out. Otherwise the full Hartle-Hawking vacuum, which is a pure 
state, should be considered. 
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in perturbation theory, an appropriate approximation scheme for the unperturbed Green's functions is also 
required. York used the Page approximation [57] for (!),„) in the Schwarzschild metric. The additional 
complication here is that while to obtain {T^ v ) the knowledge of only the thermal Fcynman Green's function 
is required, to calculate the CTP effective action one needs the knowledge of the full matrix propagator, which 
involves the Feynman, Schwinger and Wightman functions. We can put aside the technical complexity in 
the calculation of the full thermal matrix propagator G^ b [h^ v ] as our main interest is to identify and analyze 
the noise term which is the new ingredient in the backreaction problem. We have mentioned that the noise 
term gives a stochastic contribution to the Einstein Langevin equation (II. 2). We have also stated that 
this term is related to the variance of fluctuations in T^ Ul i.e., schematically, to (T^ u ). Since the Influence 
Functional or CTP formalism itself does not depend on the nature of the approximation, we will attempt to 
exhibit the general structure and project what is to be expected. 

If we denote the difference and the sum of the perturbations h^ v by [h^ v ] = h+ v — h~ v and {h^} = 
h+ v + h~ v , respectively, the influence functional form of the thermal CTP effective action may be written to 
second order in /i M „ as [56] 

T? CTP [h%] ^^rJd A x d 4 x' [h, v }(x)L^(x,x'){h^}(x') + \jd i x [h^](x)T^(x) 

+1 J d A x dV [h^](x)H^(x,x'){h a0 }(x') ~\J^x d 4 x> [h^](x)D^{x,x'){h a p}{x') 
+ l - J d A x rfV lh^}(x)N^( x ,x')lh al3 }(x'). (111.13) 

The first term is the Einstein-Hilbert action to second order in the perturbation h^ v (x) and L^ a/3 (x) is 

a symmetric local kernel, i.e. L^ a ^(x, x 1 ) — L^ a ^(x\ x). The second is a local term linear in h^ u (x). 

T^(x) represents the zeroth order contribution to (T^ v (x)} and far away from the horizon it takes the form 

of the stress tensor of massless scalar particles at temperature /3 _1 . The third and fourth terms constitute 
the remaining quadratic component of the real part of the effective action. The kernels H^' a P(x,x') and 
D^ l/ ' a/3 (x, x') are respectively even and odd in x, x' . The last term gives the imaginary part of the effective 
action and the kernel N^ v,a P(x,x') is symmetric. This is the general structure of the CTP effective action 
arising from the calculation of the traces in equation (III. 12). Of course, to write down explicit expressions 
for the non-local kernels one requires the input of the explicit form of G^ b [h^ u ] , which we have not used. In 
spite of this limitation we can make some interesting observations from this effective action. Connecting this 
thermal CTP effective action to the influence functional via equation (1.4) we see that the nonlocal imaginary 
term containing the kernel N^ ,a "(x,x') is responsible for the generation of the stochastic noise term in the 
Einstein-Langevin equation and the real non-local term containing kernel D fiv ' a,3 (x, x') is responsible for the 
non-local dissipation term. The Einstein-Langevin equation can be generated from equation (1.7) by first 
constructing the improved semiclassical effective action in accordance with (1.6) and deriving the equation of 
motion (1.7) by taking a functional derivative of the above effective action with respect to [h^] and equating 
it to zero. With the identification of the noise and dissipation kernels, one can use them to write down a 
Fluctuation-Dissipation relation (FDR) analogous to (1.9) in the context of black holes. 

IV. STRUCTURE FORMATION IN INFLATIONARY UNIVERSES 

Cosmological structure formation is a key problem in modern cosmology [59,58] and inflation [42-44] offers 
a natural solution to this problem. If an inflationary period is present, the initial seeds for the generation 
of the primordial inhomogeneitics that lead to the large scale structure observed in the present universe 
have their source in the quantum fluctuations of the inflaton field, the field which is generally responsible 
for driving inflation. Stochastic gravity provides a sound and natural formalism for the derivation of the 
cosmological perturbations generated during inflation. 
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In Ref. [16] it was shown that the correlation functions that follow from the Einstein-Langevin equation 
which emerges in the framework of stochastic gravity coincide with that obtained with the usual quantiza- 
tion procedures [60] when both the metric perturbations and the inflaton fluctuations are both linearized. 
Stochastic gravity, however, can naturally deal with the fluctuations of the inflaton field even beyond the 
linear approximation. 

Here we will illustrate the equivalence with the usual formalism, based on the quantization of the linear 
cosmological and inflaton perturbations, with one of the simplest chaotic inflationary models in which the 
background spacetime is a quasi de Sitter universe [15,16]. 

In this chaotic inflationary model [42,44] the inflaton field 4> of mass m is described by the following 
Lagrangian density 

£(<A) = lg ab v a 4>v b 4> + \m 2 <\> 2 . (iv.i) 

The conditions for the existence of an inflationary period, which is characterized by an accelerated cosmolog- 
ical expansion, is that the value of the field over a region with the typical size of the Hubble radius is higher 
than the Planck mass mp. In order to solve the cosmological horizon and flatness problem more than 60 
c- folds of expansion are needed; to achieve this the scalar field should begin with a value higher than 3mp. 
The inflaton mass is small: as we will see, the large scale anisotropies measured in the cosmic background 
radiation restrict the inflaton mass to be of the order of 10~ 6 mp. We will not discuss the naturalness of this 
inflationary model and we will simply assume that if one such region is found (inside a much larger universe) 
it will inflate to become our observable universe. 

We want to study the metric perturbations produced by the stress-energy tensor fluctuations of the 
inflaton field on the homogeneous background of a flat Friedmann-Robertson- Walker model, described by 
the cosmological scale factor a(?7), where r\ is the conformal time, which is driven by the homogeneous 
inflaton field <fi(r]) = ((f)). Thus we write the inflaton field in the following form: <f> = <f){rf) + <p(x), where 
ip(x) corresponds to a free massive quantum scalar field with zero expectation value on the homogeneous 
background metric: (ip) = 0. We will restrict ourselves to scalar-type metric perturbations because these 
are the ones that couple to the inflaton fluctuations in the linear theory. We note that this is not so if we 
were to consider inflaton fluctuations beyond the linear approximation, then tensorial and vectorial metric 
perturbations would also be driven. The perturbed metric g a b = 9ab + h ab can be written in the longitudinal 
gauge as, 

ds 2 = a 2 0y)[-(l + 2<P(x))dri 2 + (1 - 2^(x))5 ij dx i dx j ], (IV.2) 
where the scalar metric perturbations $(2;) and ^(x) correspond to Bardeen's gauge invariant variables [61]. 

A. Einstein-Langevin equation for scalar metric perturbations 

The Einstein-Langevin equation is gauge invariant, thus we can work in a desired gauge and then extract 
the gauge invariant quantities. It is given by 

G { ab ] toG(f$>) + G§(h) 8nG(f^(h)) = SttGU, (IV.3) 

where the two first terms cancel, that is G^P — 8ttG(T^) — 0, as the background metric satisfies the 
semiclassical Einstein equations. Here the subscripts (0) and (1) refer to functions in the background metric 
g a b and linear in the metric perturbation h a b, respectively. The stress tensor operator T a b for the minimally 
coupled inflaton field in the perturbed metric is: 

f ab = V a 4>V b 4> + ig ah (V c 0V c + m 2 4> 2 ). (IV.4) 

Using the decomposition of the scalar field into its homogeneous and inhomogeneous part and the metric 
gab into its homogeneous background g a b and its perturbation h a b, the rcnormalized expectation value for 
the stress-energy tensor operator can be written as 
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(IV.5) 



where the subindices indicate the degree of dependence on the homogeneous field cf> and its perturbation 
cp. The first term in this equation depends only on the homogeneous field and it is given by the classical 
expression. The second term is proportional to (<p[g]) which is not zero because the field dynamics is 
considered on the perturbed spacetime, i.e., this term includes the coupling of the field with h a b and may 
be obtained from the expectation value of the linearized Klein-Gordon equation, (o g+h — to 2 ) (p = 0. The 
last term in Eq. (IV.5) corresponds to the expectation value to the stress tensor for a free scalar field on the 
spacetime of the perturbed metric. 

After using the previous decomposition, the noise kernel N a b c d[g;x,y) can be written as 

{{tab[g;x),t cd [g;y)}) = ({i ab [g; x), i cd [g; y)})^) 2 + ({i ab [g; x), i cd [g; y)}) (w)2 , (IV.6) 

where we have used the fact that ((f) = = (>p<p<p) for Gaussian states on the background geometry. We 
consider the vacuum state to be the Euclidean vacuum which is preferred in the de Sitter background, and 
this state is Gaussian. In the above equation the first term is quadratic in (p whereas the second one is 
quartic, both contributions to the noise kernel are separately conserved since both an d & satisfy the 
Klein-Gordon field equations on the background spacetime. Consequently, the two terms can be considered 
separately. On the other hand if one treats <p as a small perturbation the second term in (IV.6) is of 
lower order than the first and may be consistently neglected, this corresponds to neglecting the last term of 
Eq. (IV.5). The stress tensor fluctuations due to a term of that kind were considered in Ref. [15]. 

We can now write down the Einstein-Langevin equations (IV. 3) to linear order in the inflaton fluctuations. 
It is easy to check [16] that the space-space components coming from the stress tensor expectation value terms 
and the stochastic tensor are diagonal, i.e. (T^) = = £y for i ^ j. This, in turn, implies that the two 
functions characterizing the scalar metric perturbations are equal: $ = * in agreement with Ref. [60]. The 
equation for $ can be obtained from the Oi-component of the Einstein-Langevin equation, which in Fourier 
space reads 

2iki(H^ k + *' fc ) = 8nG(Z 0i ) k , (IV.7) 

where ki is the comoving momentum component associated to the comoving coordinate x l , and we have 
used the definition &k(i]) = / d 3 xcxp(— ik ■ x)${r],x). Here primes denote derivatives with respect to the 
conformal time r\ and TL — a' /a. A nonlocal term of dissipative character which comes from the second term 
in Eq. (IV.5) should also appear on the left hand side of Eq. (IV.7), but we have neglected it to simplify the 
forthcoming expressions 6 . Note, however, that the equivalence of the stochastic approach to linear order in 
ip and the usual linear cosmological perturbations approach is independent of that approximation [16]. To 
solve Eq. (IV.7), whose left-hand side comes from the linearized Einstein tensor for the perturbed metric 
[60] , we need the retarded propagator for the gravitational potential , 

hrrip \ a(rj) J 

where / is a homogeneous solution of Eq. (IV.7) related to the initial conditions chosen and m 2 P = 1/G. 
For instance, if we take f(r],rf') = —d(r/ — rf)a{rf) / 'a{rj) the solution would correspond to "turning on" the 
stochastic source at 770- With the solution of the Einstein-Langevin equation (IV.7) for the scalar metric 
perturbations we are in a position to compute the two-point correlation functions for these perturbations. 



6 Such a term, which leads to an integrodifferential Einstein-Langevin equation, might not be negligible in some 
situations. Despite the apparent difficulty of dealing with an integrodifferntial equation, in this case the Einstein- 
Langevin equation can be actually transformed, after suitable manipulation, into an ordinary differential equation, 
so that the inclusion of the nonlocal term is still tractable [16] 
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B. Correlation functions for scalar metric perturbations 

The two-point correlation function for the scalar metric perturbations induced by the inflaton fluctuations 
is thus given by 

($k(v)$k>(r)'))s = (2K) 2 5(k + k')x J" d V i J" dr)2G k (mVi)G k >(r)',V2)((Zoih(m)(Zoi)k>(r)2))s- (IV.9) 

Here the two-point correlation function for the stochastic source, which is connected to the stress-energy 
tensor fluctuations through the noise kernel is given by, 

((Zoih(mmoi)- k (m))s = \{{{ki)k(m),(ioi)-k(m)})<t> v = \hhci>\mW{m)G ( k\m,m), (iv.io) 

where G^ (j/i, 772) = ({<^fe(??i), 0-k(f]2)}) is the fc-mode Hadamard function for a free minimally coupled 
scalar field in the appropriate vacuum state on the Friedmann-Robertson- Walker background. 

In practice, to make the explicit computation of the Hadamard function we will assume that the field state 
is in the Euclidean vacuum and the background spacetime is de Sitter. Furthermore, we will compute the 
Hadamard function for a massless field, and will make a perturbative expansion in terms of the dimensionlcss 
parameter m/mp. Thus we consider 

G k 1} (m,V2) = (0\{y k (vi), V -MM =2TZ (u k ( m )ul( m )) , 

with y k (r)) = a(r])ip k (r]) = a k u k (r]) + a_ k u*_ k (r]) and where u k = (2k)~ 1 / 2 e lkr '(l — i/rj) arc the positive 
frequency fc-mode for a massless minimally coupled scalar field on a de Sitter background, which define the 
Euclidean vacuum state: a k \0) = [31]. 

The assumption of a massless field for the computation of the Hadamard function is made because massless 
modes in de Sitter are much simpler to deal with than massive modes. We can see that this is, however, 
a reasonable approximation as follows. For a given mode the m — approximation is reasonable when its 
wavelength A is shorter that the Compton wavelength, A c = 1/m. In our case we have a very small mass 
m and the horizon size if -1 , where H is the Hubble constant H = a/a (here a(t) with t the physical time 
dt = adrj) satisfies that H^ 1 < A c . Thus, for modes inside the horizon A < A c and m = is a reasonable 
approximation. Outside the horizon massive modes decay in amplitude as ~ exp(— ra 2 t/ H) whereas massless 
modes remain constant, thus when modes leave the horizon the approximation will eventually break down. 
However, we only need to ensure that the approximation is still valid after 60 e- folds, i.e. Ht <~ 60, but this 
is the case since 60 m 2 < H 2 given that m ~ 10~ 6 top, and m <C H as in most inflationary models [43,59]. 

The background geometry is not exactly that of de Sitter spacetime, for which a(ij) = — (Hr/^ 1 with 
— 00 < 77 < 0. One can expand in terms of the "slow-roll" parameters and assume that to first order 
4>{t) ~ m 2 p (m/mp), where t is the physical time. The correlation function for the metric perturbation (IV.9) 
can then be easily computed; see Ref. [15,16] for details. The final result, however, is very weakly dependent 
on the initial conditions as one may understand from the fact that the accelerated expansion of de quasi-de 
Sitter spacetime during inflation erases the information about the initial conditions. Thus one may take the 
initial time to be 770 = —00 and obtain to lowest order in m/mp the expression 

(*fc(»7)*k' (»/)>- - S77 2 (m/mp) 2 k- 3 (2ir) 3 8(k + k') cosk^ - r/). (IV.ll) 

From this result two main conclusions are derived. First, the prediction of an almost Harrison-Zel'dovich 
scale-invariant spectrum for large scales, i.e. small values of k. Second, since the correlation function is of 
order of (m/mp) a severe bound to the mass m is imposed by the gravitational fluctuations derived from 
the small values of the Cosmic Microwave Background (CMB) anisotropics detected by COBE [60]. This 
bound is of the order of (m/mp) <~ 10~ 6 . 

We should now comment on some differences with those works in Ref. [63] which used a self-interacting 
scalar field or a scalar field interacting nonlinear ly with other fields. In those works an important relaxation 
of the ratio m/m p was found. The long wavelength modes of the inflaton field were regarded as an open 
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system in an environment made out of the shorter wavelength modes. Then, Langevin type equations were 
used to compute the correlations of the long wavelength modes driven by the fluctuations of the shorter 
wavelength modes. In order to get a significant relaxation on the above ratio, however, one had to assume 
that the correlations of the free long wavelength modes, which correspond to the dispersion of the system 
initial state, had to be very small. Otherwise they dominate by several orders of magnitude those fluctu- 
ations that come from the noise of the environment. This would require a great amount of fine-tuning for 
the initial quantum state of each mode [16]. We should remark that in the model discussed here there is no 
environment for the inflaton fluctuations. The inflaton fluctuations, however, are responsible for the noise 
that induce the metric perturbations. 
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